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Abstract 

Let (L; n, U) be a finite lattice and let n be a positive integer. A 
function / : L" ^ R is said to be submodular if /(a fl 6) + f{a U 6) < 
f{a) + f{b) for all a,b £ L". In this paper we study submodular functions 
when L is a diamond. Given oracle access to / we are interested in finding 
X £ L" such that f(x) — miUj^gLn f{y) as efficiently as possible. We 
establish 

• a min-max theorem, which states that the minimum of the submod- 
ular function is equal to the maximum of a certain function defined 
over a certain polyhedron; and 

• a good characterisation of the minimisation problem, i.e., we show 
that given an oracle for computing a submodular / : L" — > Z and an 
integer m such that miuajgin fi^) = "^-i there is a proof of this fact 
which can be verified in time polynomial in n and maxt^L^ log 

and 

• a pseudo-polynomial time algorithm for the minimisation problem, 
i.e., given an oracle for computing a submodular f : L" Z one 
can find miutgL" fit) in time bounded by a polynomial in n and 



maxtgL" 

1 Introduction 



Let y be a finite set and let / be a function from 2^ to M. The function 
/ is said to be submodular if f{A U B) + f{A D B) < f{A) + f{B) for all 
B <Z V . In the sequel we will call such functions submodular set Junctions. 
Submodular set functions shows up in various fields including combinatorial 
optimisation, graph theory [7], game theory |30| . information theory ^13j and 
statistical physics [T] . Examples include the cut function of graphs and the rank 
function of matroids. There is also a connection between submodular function 
minimisation and convex optimisation. In particular, submodularity can be seen 
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as a discrete analog of convexity [51 [23] . We refer the reader to [HI [HI [M] for a 
general background on submodular set functions. 

Given a submodular set function / : 2^ ^ M there are several algorithms 
for finding minimisers of /, i.e., finding a subset X C V such that f{X) = 
minycy f(Y), in time polynomial in \V\. The first algorithm for finding such 
minimisers in polynomial time is due to Grotschel et al. [IT]. However, this al- 
gorithm is based on the Ellipsoid algorithm and hence its usefulness in practise 
is limited. Almost two decades later two combinatorial algorithms were found 
independently by Schrijver [55] and Iwata et al. [IH] • More recently the running 
times have been improved. The currently fastest strongly polynomial time al- 
gorithm is due to Orlin [5S] and the fastest weakly polynomial time algorithm 
is due to Iwata [17]. In these algorithms the submodular set function is given 
by a value-giving oracle for / (i.e., presented with a subset X C V the oracle 
computes f{X)). 

In this paper we investigate a more general notion of submodularity. Recall 
that a lattice is a partially ordered set in which each pair of elements have a 
least upper bound (join, U) and a greatest lower bound (meet, □). Given a 
finite lattice C (all lattices in this paper are finite) and a positive integer n we 
can construct the product lattice Meet and join for £" are then defined 
coordinate-wise by meet and join in £. We say that a function ft. : £" — > M is 
submodular if h{a nb) + h{a U b) < h{a) + hib) for all a, b e Note that 
the subsets of V can be seen as a lattice with union as join and intersection as 
meet (this lattice is a product of the two element lattice) . Hence, this notion 
of submodularity is a generalisation of submodular set functions. For a fixed 
finite lattice C we are interested in the submodular function minimisation (SFM) 
problem: 

Instance: An integer n>l and a submodular function / on 
Goal: Find a; G such that f{x) — vaiTiy^c^ fiv)- 

Following [22] we denote this problem by SFM(/:). SFM(Z:) is said to be 
oracle-tractable if the problem can be solved in time polynomial in n (provided 
that we have access to a value-giving oracle for / and that we can assume that 
/ is submodular, i.e., it is a promise problem). This definition naturally leads 
to the following question: is SFM(£) oracle-tractable for all finite lattices £? 
(This question was, as far as we know, first asked by Cohen et al. [3].) 

Schrijver [28] showed that given a sublattice S of 2^ (i.e., S* C 2^ and 
for any X,Y G S we have X 11 Y, X U Y G S) and submodular function / : 
5 — > R a minimiser of / can be found in time polynomial in n. In particular, 
this implies that for any distributive lattice £ the problem SFM(£) is oracle- 
tractable. Krokhin and Larose [22] showed that certain constructions on lattices 
preserve oracle-tractability of SFM. In particular, they showed that if AT is a 
class of lattices such that SFM(£) is oracle-tractable for every C e X, then so 
is SFM(/^') where C is a homomorphic image of some lattice in AT, a direct 
product of some lattices in X , or contained in the Mal'tsev product X o X . We 
will not define these constructions here and refer the reader to [22] instead. 

A lattice C is a diamond if the elements of the lattice form a disjoint union 
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Figure 1: The five element diamond. 

of {0£, l^} and A, for some finite set A such that \A\ > 3. Here 0£ is the 
bottom element of C, and Ic is the top element of C, and all elements in A 
(called the atoms) are incomparable to each other. See Figure [1] for a diagram 
of the five element diamond. We want to emphasise that diamonds have a 
different structure compared to the lattices defined by union and intersection. In 
particular, diamonds are not distributive, that is they do not satisfy xn{yUz) = 
(x ny) Li {x n z) for all x,y,z £ C. We will denote the diamond with k atoms 
by Mk- In Sections [H [5] and [H] the complexity of SFM{Mk) is investigated. In 
the approach taken in this paper the difficult case is fc = 3 — the proofs for the 
A; = 3 case generalises straightforwardly to an arbitrary k. We note that none 
of the diamonds are captured by the combination of the results found in [22l EH] 
(a proof of this fact can be found in [22]). 

Results and techniques. The first main result in this paper is a min-max 
theorem for SFM(7Wfc) which is stated as Theorem 14.31 This result looks quite 
similar to Edmonds' min-max theorem for submodular set functions 'F' (we 
present Edmonds' result in Section [2]). The key step in the proof of this result 
is the definition of a certain polyhedron, which depends on /. 

The second main result is a good characterisation of SFyi{A4k) (Theo- 
rem [5^. That is, we prove that given a submodular f : Z, and integer m 
such that mina;g£n f{x) = m, there is a proof of this fact which can be verified 
in time polynomial in n and maxyg£n log|/(y)| (under the assumption that / 
is submodular). This can be seen as placing SFM(7Wfc) in the appropriately 
modified variant of NP n coNP (the differences from our setting to an ordinary 
optimisation problem is that we are given oracle access to the function to be 
minimised and we assume that the given function is submodular). The proof of 
this result makes use of Caratheodory's theorem and of the known polynomial- 
time algorithms for minimising submodular set functions. We also need our 
min-max theorem. 

The third result is a pseudo-polynomial time algorithm for SFM(A^fc) (see 
Section[6]). We show that SFM(A^fc) can be solved in time polynomial in n and 
maxtgAij l/(*)l- The main part of the algorithm consists of a nested applica- 
tion of the Ellipsoid algorithm. We also need to prove that the polyhedrons we 
associate with submodular functions are 1/2-integral. An interesting and chal- 
lenging open problem is to construct an algorithm with running time polynomial 
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Figure 2: A lattice which can be shown to admit a pseudo-polynomial time 
algorithm for the submodular function minimisation problem. This lattice is 
a Mal'tsev product of a diamond and the two element lattice. By the results 
in this paper and the constructions in [22j this lattice gives a new tractable 
constraint language for Max CSP. 

in n and maxt^Mi log |/(*)|- 



Our results applies to diamonds, however, as mentioned above, in [22| two 

constructions on lattices (Mal'tsev products and homomorphic images) are shown 
to preserve tractability results for SFM. By combining these constructions with 
the results in this paper onegets tractability results for a much larger class of 
lattices than just diamonds. □ In particular, by the results in this paper there is 
a pseudo-polynomial time algorithm for minimising submodular functions over 
products of the lattice in Figure [21 

Connections to other problems. Minimising submodular functions on cer- 
tain modular non-distributive (the diamonds are modular and non-distributive) 
lattices has applications to canonical forms of partitioned matrices [TSl 1^ . 
Finding a polynomial time algorithm for minimising submodular functions on 
these lattices was mentioned as an open problem in [1_9J . 

The notion of submodular functions over arbitrary finite lattices plays an 
important role in the complexity of the maximum constraint satisfaction prob- 
lem (Max CSP). This connection was first observed in f3] and in later papers 
the connection was explored further [5] |21| . The connection between submodu- 
lar function minimisation and Max CSP is that by proving oracle-tractability 
for new lattices for the SFM problem implies tractability results (solvability in 
polynomial time) for certain restrictions (so called constraint language restric- 
tions) of Max CSP. By constructing algorithms for SFM with running times 
bounded by a polynomial in n and maxteM]^ we do in Section [51 one 

gets solvability in polynomial time for the unweighted variant of Max CSP 
(with appropriate restrictions). Providing good characterisations of SFM(£), 
as we do in Section [5l implies coNP containment results for Max CSP (with 
appropriate restriction). As Max CSP is trivially in NP we get containment 
in NP n coNP for these restrictions. We refer the reader to [31 [^ for further 
details regarding the connection between SFM and Max CSP. 

^In 1221 these constructions are shown to preserve oracle-tractability and not solvability in 
pseudo-polynomial time. However, it is straightforward to adapt the proofs to the pseudo- 
polynomial case. 
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In [22] it is shown that the restrictions of Max CSP which one gets from 
the diamonds can be solved in polynomial time. This means that the results 
for the diamonds in this paper does not directly imply new tractability results 
for Max CSP. However, as mentioned in the previous section one can combine 
the results in this paper with the lattice constructions in [55] to get tractability 
results for a larger class of lattices which implies tractability results for new 
constraint language restrictions of Max CSP. (We again refer to Figure [2] for 
an example of such a lattice.) 

There is also a connection between SFM over lattices to the Valued Con- 
straint Satisfaction Problem (VCSP). See [1] for more information on VCSP. The 
connection is very similar to the connection to Max CSP, proving tractabil- 
ity results for new lattices for SFM implies new tractable restrictions of VCSP. 
For VCSP there was, before the results in this paper, no known non-trivial 
algorithms for the restrictions one obtains from the diamonds. 

We note that Raghavendra [SS] recently proved almost optimal results for 
the approximability of MAX CSP for constraint language restrictions, assuming 
that the unique games conjecture (UGC) holds. However, for the cases which are 
solvable to optimality the results in [26] gives us polynomial-time approximation 
schemes. This should be compared to the connection described above which 
gives polynomial time algorithms for some of these cases. 

Organisation. This paper is organised as follows, in Section [5| we give a short 
background on submodular set functions, in Section[3|we introduce the notation 
we use, in Section [4] we prove our first main result — the min-max theorem for 
submodular functions over diamonds. The good characterisation is given in 
Section [5] In Section [6| where we give the pseudo-polynomial time algorithm 
for the minimisation problem. Finally, in Section [7| we give some conclusions 
and open problems. 

2 Background on Submodular Set Functions 

In this section we will give a short background on Edmonds' min-max theorem 
for submodular set functions. This result was first proved by Edmonds in [6 , but 
see also the surveys [18l[24]. Let y be a finite set. For a vector x e MX (i.e., x 
is a function from V into M) and a subset Y CV define x{Y) = J2yeY ^iu)- We 
write a; < if x{v) < for all v G V and x^ for the vector in which coordinate 
V has the value min{0, x{v)}. Let / be a submodular set function / : 2^ ^ M 
such that /(0) — (this is not really a restriction, given a submodular function 
g we can define a new function g'{X) = g{X) — (?(0), g' satisfies g'(0) = and is 
submodular). The submodular polyhedron and the base polyhedron defined by 

P(/) = {cc e I Vr C 1/, x(Y) < f(Y)}, and 
B{f) = {xeR''\xE P{f),x{V) - f{V)} 



5 



3 PRELIMINARIES 



often play an important role in results related to submodular set functions. 
Edmonds [6l proved the following min-max theorem 



In Section!?] we give an analog to ([T]) for submodular functions over diamonds. 

3 Preliminaries 

For a positive integer n, [n] is the set {1,2,..., n}. Given a lattice (L, PI, U) and 
x,y € L we write x Q y if and only if xFly = x (and hence xUy = y). We write 
X \Z y ii X Q y and x ^ y. As mentioned in the introduction, given a positive 
integer n, we can construct the product lattice from L. The top and bottom 
elements of L" are denoted by l^n and O^n, respectively. We write x -< y if x 
is covered by y (that is, if x n y, and there is no z G L such that x \Z z \Z y). 

Recall that the diamonds are modular lattices (the rank function p is defined 
by p{Om) — 0; — 1 for ail a G A and p{1m) = 2). As direct products of 
modular lattices also are modular lattices it follows that direct products of 
diamonds are modular lattices. 

For a set X we let rNx-^ be the set of functions mapping [n] x X into 
R. Such functions will be called vectors and can be seen as vectors indexed by 
pairs from [n] x X. For x,y d rNx^ and a G M we define ax,x + y,a;~ £ 
rNx-^ as {ax){i,x) = ax{i,x), {x + y){i,x) — x{i,x) +y{i,x), and x^{i,x) = 
min.{0,x{i,x)} for all i G [n] and x G X, respectively. If x{i,x) < for all 
i G [n] and x X we write a; < 0. For i G [n] we use x{i) to denote the function 
x' G M.-^ such that a;(i,a;) = x'{x) for all a: G X. 

For i G [n] and a G A let Xi,a G rI"!^"^ be the vector such that Xi,a{h a) — 1 
and Xi,a{i\a') = for {i',a') ^ (»,«)• (So Xj^a is the unit vector for the 
coordinate (i,a).) Similarly, we use Xi to denote the vector X]aeA/ti,a- For a 
vector X G RI"!^'^ and tuple y G A^" we define 



(This should be compared to how applying a vector to a subset is defined for 
submodular set functions, see [3].) For x, x' G RWx-4 denote the usual scalar 
product by {x,x'), so 



nun /(X) = max{cc(\/) | x G P{f),x < 0} 
= max{a;-(V") | x G B{f)}. 



(1) 



n 




n 
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Let / be a submodular function on A^" such that /(Om") > 0. We define 
PmU) and BM{f) as foUows, 



Due to the definition of g it is not hard to see that PmU) is a polyhedron. Note 
that if t contains at least one then t induce more than one linear inequality. 
If t contains no 1m i then t only induce one linear inequality. In general, a tuple 



denote the set of all vectors e G RI"!^'^ such that e represents an inequality 
induced by t (that is, an inequality of the form (e,cc) < /(t), where e G 
Given a vector x G PAi{f) we say that a tuple t G A^" such that x{t) = f{t) is 
X -tight. 

We will also need the following definition. 

Definition 3.1 (Unified Vector for Diamonds). A vector x G K"^ is unified if 
there is an atom p ^ A such that 

• if x,y A,x,y ^ p, then x{x) = x{y); and 

• if a ^ A, then x{p) > x{a). 

We extend the definition of unified vectors to the vectors in RWX'^ f^y guying 
that X G RI"!^"^ is unified if x ^ x{i, x) is unified for each i G [n]. 

If the submodular inequality is strict for all incomparable pair of elements 
then we say that the function is strictly submodular. 

4 A Min-Max Theorem 

The main results in this section are Theorem l4.3l and Theorem l4.5l We start by 
a lemma which shows that BM{f) is non-empty for any submodular function 
which maps the bottom of the lattice to a non-negative value. 

Lemma 4.1. Let f : yW" ^ R 6e submodular such that f{0) > 0. There is a 
vector X G RI"!^-^ such that 

• X is unified; and 

• x{vi) ~ f{vi) for all i G [n\; and 

• x{vi[i + 1 = Pi+i]) = f{vi[i + l = Pi+i]) for alii G {0, 1, . . . , n — 1}, where 
for i G [n], pi is the atom in Definition \3.1\ for the vector x i-^ x{i,x). 

Furthermore, if f is integer-valued, then x can be chosen to be integer-valued. 




with m occurrences of 1m induces 
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Proof. Given a submodular / : yVf" R we will construct a vector x which 
satisfies the requirements in the lemma. To do this we define a sequence of 
atoms Pi for i G [n] inductively. To start the inductive definition let pi G 
max arg^g^ fi'^oi^ = a]) and set x{l,pi) — f{vo[l ~ pi]). For the general case, 
choose Pi € A so that 

Pi G max arg/(t!i[i + 1 a]). 

aeA 

For i G [n] set 

x{i + l,p^+l) = f{vi[i + 1 = p,,+i]) - f{vi), (2) 
and for a G A, a 7^ Pi+i set 

x{i + 1, a) - /(-yi+i) - f{vi\i + 1 = (3) 

Claim A. If a G A, then x{i + > x{i + 1, a). 

Assume, without loss of generality, that a ^ Pi+i- We now get 

f{vi[i + l=pi+i]) + f{vi[i + l^a]) < 
2f{vi[i + 1 = Pi+i]) 

where the first inequality holds due to the submodularity of / and the second 
inequality follows from our choice of Pi+i- This is equivalent to 

/K+i) - fivi[i + 1 = Pr+l]) < f{Vi[i + 1 = - f{Vi) 

which is what we wanted to prove. □ 
For i G [n] and j G {1, 2} we define Cij as Ci^i = pi and Ci^2 = 1m- n 
Claim B. x{vi[i + 1 = q j]) = f{vi[i + \ = a j]) for all (i, j) G {0, 1, . . . , n — 
l}x{l,2}. 

We prove this by induction over the pairs (i, j) ordered lexicographically (so 
< (^'ij') if and only li i < i' or (? = i' and j < j')). With the pair 
we associate the tuple Vi[i + 1 = Cij]- Note that < if and only if 

+ 1 = Ci+ij] C + 1 = Ci/+ij/]. As pi G max arg„g^/(vo[l = a]) the 
claim clearly holds for ~ (0,1)- Now assume that it holds for all pairs 

such that < If J = 1 then the next pair is (i,2) and we get 

x{vi[i + 1 = Q.2]) = x{vi[i + 1 = Pi+i]) + x{i + l,a) 

= f{vi[i + 1 = Pi+i]) + f{vi+i) - f{vi[i + 1 = p,+i]) 

= f{Vi+l)- 

Here the first inequality follows from the definition of x{-) and Claim A. The 
second equality follows from the induction hypothesis and ([3]) . If j = 2 the next 
pair is (i + 1, 1) and we get 

x{vi+i[i + 2 = Ci+2,1]) = x{vi+i) + x{i + 2,pi+2) 

= f{vi+i) + fivi+i[i + 2=p,+2]) - fivi+i) 
= fivi+i[i + 2=^+2]) 
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As above the first equality follows from the definition of x{-) and Claim A. The 
second equality follows from the induction hypothesis and ([3]). □ 
By Claim A it follows that x is unified. By Claim B x satisfies the second 
condition in the statement of the lemma. It is easy to see that if / is integer- 
valued, then so is a;. □ 

Lemma 4.2. Let f : A^" ^ K. fee suhmodular such that /(O) > 0. Let x he a 

vector in m1"1'*"^. Lf for each i G [n] there is an atom pi such that 

• for all i G [n] we have x(vi) ~ f(vi), and 

• for all i G {0, 1, . . . , 71—1} we have x{vi[i+l = Pi+i]) = f{vi[i+l = Pi+i]), 
then X G Pm(/)- 

Proof. For i G [n] and j G {1, 2} we define Cij as follows c,;.i = pi and Ci^2 = ^m- 
We will prove by induction that x{y) < f{y) for all y G As in the proof 

of Claim B in Lemma |4J] the induction will be over the pairs {0, 1, . . . , rt — 1} x 
{1,2} ordered lexicographically. With the pair we associate the tuples y 
such that y ^ Vi[i + 1 = Pij]- 
As 

x{vo) - xiOMr.) ^ and /(0a4") > 

and 

x{vo[l ^pi]) ^ /(fo[l = Pi]) 

the statement holds for the pair (0, 1) (which corresponds to y □ [1 = pi]). 
Let i G {0, 1, . . . , n - 1}, j G {1, 2}, and y G y C + 1 = c,:+ij] and 
assume that the inequality holds for all y' G Al" such that y' C Vii [j' + 1 = 
Ci'+i_j'] where {i' is the predecessor to the pair {i,j). We will prove that the 
inequality holds for b1\ y ^ Vi[i + 1 — c^+ij]. 

To simplify the notation a bit we let fc = i + 1 and y = y{k). If y = O^vi we 
are already done, so assume that y ^ O^. If y = Pk let c = O^vi, if y A, y ^ pk 
let c = Pk and otherwise, if y = Im let c = pk- Now, 

x{y) < x{vi[k 

< x{vi[k 

< f{Vi[k 

< f{y)- 

The first inequality follows from the supermodularity of x. The second in- 
equality follows from the induction hypothesis and the fact that y Fl c n y and 
y n c G {0;\4,pfc}. The third inequality follows from y U c, c G {OM^PkAM} 
and the assumptions in the statement of the lemma. Finally, the last inequality 
follows from the submodularity of /. □ 
In the proof of Lemma HTTl the vector x G K^"! ^'^ is constructed with a greedy 
approach — we order the coordinates of the vector, [n] x A, in a certain way and 
then set each component to its maximum value subject to the constraints given 



= y U c]) — x{vi[k = c]) + x{y[k ~ yH c]) 
= y U c]) - x{vi[k = c]) + f{y[k = y n c]) 
= y U c]) - f{vi[k - c]) + f{y[k - y n c]) 
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in the definition of Bi\i (/) . The greedy algorithm does not solve the optimisation 
problem for PmU)- As an example, let Mz = ({Otvi, l^vi, a, 6, c}, U) be a 
diamond and let / : Mz — > M be defined as /(Ox) — 0, /(a) = f{h) = 
/(c) = /(l^) = 1. The function / is submodular. Now let c € IRI^]^^ and 
c(l, a) = c(l, 6) = c(l, c) = 1. From the greedy algorithm we will get a vector 
X S RI-^I^^ such that a;(l, a) = 1 and a;(l, 6) = a;(l, c) = (or some permutation 
of this vector). However, the solution to max(c, y),y G -Pm(/) is y(l,a) = 
y(l,6) = y(l,c) = 1/2 and 3/2 = (c, y) > (c, a;) = 1. This example also shows 
that the vertices of PM{f) are not necessarily integer valued. This should be 
compared to submodular set functions, where the corresponding optimisation 
problem is solved by the greedy algorithm. [24] 

Given an algorithm which solves the optimisation problem over Pm{J) in 
time polynomial in n we can use the equivalence of optimisation and separation 
given by the Ellipsoid algorithm to solve the separation problem for PmH) in 
polynomial time. With such an algorithm we can decide if g -Pm(/) or not 
and by a binary search we can find a minimiser of / in polynomial time. So a 
polynomial time algorithm for the optimisation problem over PmH) would be 
desirable. (The approach outlined above can be used to minimise submodular 
set functions, see [IT] or, e.g., [T^].) We present a pseudo-polynomial algorithm 
for the optimisation problem in Section [H] which uses this technique. 

We are now ready to state the two main theorems of this section. 

Theorem 4.3. Let / : -^M. he a submodular Junction such that /{Om") = 
0, then 

min f{x) = max{z{lM.rz) I z G PMif),z < 0, z is unified } . 

More over, if f is integer-valued then there is an integer-valued vector z which 
maximises the right hand side. 

Proof, li ze PmU) and 2 < then 

z{1m-) < z{y) < f{y) 

for any y G M'\ Hence, LHS > RHS holds. Consider the function /' : M" M 
defined by 

f'{x) = niin/(y). 

y\Zx 

Then Pm(/') C Pm(/)- 
Claim A. /' is submodular. 

Let x',y' G Ai" and let x C x',y □ y' be tuples such that f'{x') ~ f{x) 
and /'(y') = /(y). Now, 

fix') + f'{y') = fix) + f{y) > fix n y) + f{x U y) > f'{x' H y') + f'{x' U y') 

where the first equality follows from the definition of /', x and y, the first 
inequality follows from the submodularity of / and the second inequality from 
the definition of /' and xHy \—x' \ly' and xUy \—x' Uy' . □ 
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Claim B. For any z € Puif) we have z <0. 

As /(Oa<>.) = we have f'{x) < for any x e X". For i G [n] and a e A 
define U^a G A^" such that *i,a(j) = Om for j E [n],j ^ i and ti,a(«) a- It 
follows from 2; £ -Pa/(/') that we have z{ti^a) — z{i,a) < f'{ti^a) < for any 
a € A and i G [n]. □ 
Claim C. Any z e Sm(/') C Pm(/') satisfies ^(1^.) = /'(Ia^-)- 

Follows from the definition of BM{f') □ 

Finally, /'(1a<") = min^jg^n f{x) which follows from the definition of /'. 
From Lemma ini and Lemma lT^ it now follows that LHS < RHS holds. To prove 
the existence of a integer valued vector, note that the vector from Lemma 14.11 
is integer valued if /' is integer valued and /' is integer valued if / is integer 
valued. □ 

We can reformulate Theorem 14.31 to relate the minimum of a submodular 
function / to the maximum of a certain function defined over the polyhedron 
{x e Pm{/) I a; < 0}. To do this we define a function S : mI"]^^ ^ M as follows 

n 

Six)—} mina;(i, a) + maxccfz, a). 

^ — ' aGA aGA 

We then get the following corollary. 
Corollary 4.4. 

min f{y) = max {S{z) \ z e PmU). z<G]. 



Proof. Follows from Theorem 14.31 by two observations. If z is unified, then 
z(l7V4") = S{z). Furthermore, any vector z can be turned into a unified vector 
z' such that z' < z and S{z) = (To construct z' from z, for each 

i € [n], choose some pi G max arg^^g^ a) and let z'{i,pi) — z{i,pi) and for 
a <E A, a ^ Pi let z'{i, a) = miuag^ z{i, a).) □ 
One might ask if there is any reason to believe that the min~max character- 
isation given by Theorem 14.31 is the "right" way to look at this problem. That 
is, can this min-max relation give insight into the complexity of minimising 
submodular functions over diamonds? Theorem l4.3l is used in Section [S] to get a 
good characterisation of submodular function minimisation over diamonds, so it 
certainly gets us somewhere. In Section [5] we present a pseudo-polynomial time 
algorithm which uses Pnif), but it does not use Theorem 14.31 Additionally, 
Theorem 14. 31 is in some sense fairly similar to ([1]). In particular, in both cases 
the vectors are functions from the atoms of the lattices to the real numbers 
and when a vector is applied to a tuple (or a subset) it is computed as a sum 
over the coordinates of the vector and the tuple. Furthermore, in this sum the 
bottom of the lattice {Qm in the diamond case and in the set case) do not 
contribute to the sum. There are of course differences as well. The most obvi- 
ous one is, perhaps, that there is no element in the set case analogous to 1m 
in the diamond case. Considering that, as far as we know, all combinatorial 
algorithms for submodular set function minimisation is based on ^ and the 
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similarity between Theorem 14.31 and H]) one could hope that Theorem 14 . 31 could 
be the basis for a polynomial time combinatorial algorithm for SFM(AI). 

The following theorem is an analog to the second equality in Edmonds' min- 
max theorem for submodular set functions ([T]) . 

Theorem 4.5. Let f : A4" -^R be a submodular function such that /(O^vfn) = 
0, then 

min f{x) = max|2;(l^n) I z G PM{f),z < 0, z is unified] 

— max {a;^(l^Ti) | x E BM{f), x is unifiedj . 

Proof. We prove that 

max {2;(ljv4Ti) | z G PM{f),z < 0,2; is unified} = 
max{a;^(l^n) | x G B]\f{f),x is unified} . 

The result then follows from Theorem l4.3l 

Let a; be a vector which maximises the right hand side. It is clear that 
x~ G PmH), < 0, and that x~ is unified. It follows that LHS > RHS. 

Conversely, let 2; be a vector which maximises the left hand side. We will 
define a sequence of vectors Xq,Xi, . . .. We start with Xq = z and for j > we 
define ajj+i from xj according to the construction below. 

1. If there is some i G [n] and p G max arg^^^ x{i, a) such that a' > where 

a' = max{Q: G M | a; + axi,p e PAiif)}, 
then let Xj^i — Xj + a' ■ Xi,p- 

2. Otherwise, if there is some i G [n] and p G max arg^^^ x{i, a) such that 
a' > where 

a' = max{a e M. \ Xj + a ■ {xi - Xi,p) e PAiif)}, 

then let a be some atom distinct from p, let m — min{a', x{i,p) —x{i, a)}, 
and let Xj+i =Xj +m- (xi - Xi,p)- 

We make four observations of this construction. 

• If we reach the second step, then max arg^g^ x{i, a) is a one element set. 

• For every j the vector xj is unified. 

• For every j, Xj^i > Xj. 

• For every j, xj G P{f)- 
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These observations all follows directly from the construction above. It is not 
hard to convince oneself that there is an integer m such that Xm — Xm+i (and 
thus all vectors constructed after m are equal). To see this, note that for a fixed 
i G [n] if some atom a is increased in step 1, then this atom will not be increased 
again at coordinate i. Let y denote the vector Xm- 

Note that a;j_|_i~(l7vi") > a;j~(l^>>) for all j. Hence in particular y~ (l^vi") > 
z{1m")- As we have already proved that LHS > RHS it now remains to prove 
that y e BmH)- As we already know that y £ Pa/(/) this reduces to proving 
y(W) -/(1a4"). 

Let p be a tuple such that for i G [n] we have p{i) = max^gA y(*, a). As 
y = Xm = Xm+1, it follows that for each k 6 [n] there is an atom a G A,a 
p[k) and tuples tfc,<^ G M'^,p{k) Q tk{k),a C t^(fc) such that tk and are 
y-tight. Now let, 

ke[n] 

As y G PmH) it follows from Lemma l5.ll that y(t) = f{t). Note that for 
each k E [n] we have {t^ U i^)(fc) = Im, it follows that t = l^^n and hence 
y e BmU). We conclude that LHS < RHS. □ 

5 A Good Characterisation 

In this section we show that there are membership proofs for Pa/(/) which can 
be checked in time polynomial in n. By using Theorem 14. 31 this will lead to the 
existence of proofs that can be checked in time polynomial in n of the fact that 
a certain tuple minimises a submodular function. The following lemma states 
that if a and b arc a;-tight, then so are a Fib and a Lib. This simple result will 
be used repeatedly in the subsequent parts of the paper. 

Lemma 5.1. Let f : A^" M. be a submodular function. Let x e PmU) be a 
vector and let a,b G be x-tight tuples. Then, aUb and aHb are x-tight. 

Proof. 

x{a Ub) + x{a n 6) < /(a U b) + /(a n 6) < /(a) + f{b) = x{a) + x{b) 

The first inequality follows from x G Phiif), the second follows from the sub- 
modularity of /. The equality follows from the assumptions in the lemma. 
Note that x{a) + x{b) < x{a Ub) + x{a n b). Since x{a U b) < /(a U b) and 
x(anb) < /(an 6), it follows that x{aUb) = f{aUb) and x{anb) = /(an 6). 

n 

The following lemma is an important part of the main result in this section. 

Lemma 5.2. Let c G R["1^^ and assume that x maximises {x, c) over Puif)- 
Furthermore, assume that a,b G A^", a C 6 are x-tight and for all t G Ad" such 
that a n t \Z b the tuple t is not x-tight. Then, there is at most one coordinate 
i G [n] such that a(i) — Om cind b{i) — 1m- 
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Proof. Assume that there is another coordinate j S ["], j ^ « such that a{j) = 
Om and b{j) ~ 1m- We can assume, without loss of generahty, that 

x^A xeA 

Let 5 > Q and let x' = x + 5xi ~ ^Xj- We cannot have x' g Pm(/) for any 
6 > 0, because then x is not optimal. As x' ^ Pm(/) there is some a;-tight tuple 
t e X" such that (t(i) G A andt(j) = O^i) or {t{i) = 1m and e {Oai}U^). 
In either case, it follows from Lemma that t' = (6 n f ) U a is a;-tight, which 
is a contradiction as a IZ t' C b. □ 
The key lemma of this section is the following result. We will use this 
lemma together with Lemma 15.21 in the proof of the main result of this section 
(Theorem [511). 

Lemma 5.3. Let n be a positive integer and let f : A^" M be submodular 
which is provided to us by a value-giving oracle. Let x G M^"!^"^ and a,b Cz Al" 
such that a \—b, a is x-tight, and there are at most k coordinates i G [n] such 
that a(i) — Ox and b(i) = 1m- Under the assumption that for all t Q a we 
have x{t) < f{t) it can be verified in time O (n'^^'^) that x{y) < f{y) holds for 
all y Q b, for some fixed constant c. 

Proof. Let / C [n] be the set of coordinates such that i G / if and only if 
a{i) = Om and b{i) = 1m and let J — {j £ [n] \ a{j) ^ b{j),j ^ I}. Let 
Z = {z e A4" I Vi ^ / : z{i) = Om}- For a subset Y = {yi, y2, - - - , Vm} of J 
and z ^ Z define gz : 2'^ R as 

gziY) = f{a[yi = b{yi), ...,y^ = b{y,n)] U z). 

We claim that gz is a submodular set function. Let z E Z and let C = 
{ci, C2, . . . , Cfe} and D = {di, ^2, . . . , d;} be two arbitrary subsets of J. De- 
fine c,d £ A^" as a[ci — b(ci), . . . , Cfc = b{ck)] U z and a[di — b{di), . . . ,di = 
b{di)] U z, respectively. We now get 

gz{C)+gziD)^f{c) + f{d) > f{cnd) + f{cUd)^gz{CnD)+gz{CUD). 

Hence gz is submodular for each z E Z. For a subset Y ~ {yi, y2, ■ • • , J/m} of J 
define /i^ : 2'^ ^ R as 

hziY) ^ x{a[yi = %i), . . . , y„ = U z). 

We claim that —hz is a submodular set function for each z G Z. As above, 
let C = {ci, C2, . . . , cj.} and D — {c?i, d2, . - - , di} be two arbitrary subsets of J 
and let c — a[ci — b(ci), - - - ,Ck = b(cfc)] U z and d = a[di = b{di), - - - ,di = 
b{di)] U z, then 

hz{C)+hz{D) = x{c)+x{d) < x{cnd)+x{cUd) = hz{C n D) + hz{C U D). 
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Hence, —hz is submodular. Let Y = {yi, 1/2, ■ ■ ■ , ym} be an arbitrary subset of 
J and let z G Z. For a fixed k the inequalities 

x{a[yi = %i), ...,ym = b(2/m)] U z) < 
f{a[yi = 6(2/1), . . . , y„ = U 2;) 

0<5.(r)-/^4i") (4) 

can be verified to hold for every Y <Z J and z E Z in time polynomial in n as, 
for each x E X, the RHS of is a submodular set function in Y. Conversely, 
if ([U does not hold for some Y C J and x G X, then there is a tuple t \Z b 
such that x{t) ^ /(*). To verify that (g]) holds for all F C J and 2; e Z find 
the minimum value of the RHS of ^ for each z E Z and compare it to (note 
that \Z\ only depends on k and \A\). This can be done in time polynomial in n 
by one of the polynomial time algorithms for submodular function minimisation 
(see, e.g., [T21II11I1H] for descriptions of these algorithms). 

Let y G A^" be a tuple such that y Q b. Note that if a ^ y, then it 
follows from (|4]) that x{y) < f{y). For the sake of contradiction, assume that 
x{y) ^ /(?/)• By the submodularity of / we get 

f{aUy) + fiany)<f{a) + f{y). (5) 

AsaCaUyCbit follows from ^ that x{a Liy) < f{a U y). Furthermore, 
y n a CI a so by the assumptions in the lemma x{a n y) < /(a □ y). By the 
choice of a and y we get x{a) = f{a) and /(y) < x{y). It follows that 

x{aUy)+x{any)<f{aUy) + J{any) (6) 

and 

f{a) + f{y)<x{a)+x{y). (7) 

But 

x{a) + x{y) < x{aUy) + x{a\ly) (8) 

so we get a contradiction by combining ©, ©, ([7]), and ([H]). □ 
Before we prove the main result of this section we need a few basic facts 
about polyhedrons. Let P C M" be a polyhedron. The lineality space of P, 
denoted by lin.space P, is the set of vectors x such that there is a vector y E P 
and Aa; + y G P for all A £ R. The characteristic cone of P, denoted by 
char. cone P, is the set of vectors x 6 M" such that for all y G P and A > we 
have Ace + y e P. 

Given a submodular function /, it is not hard to see that the characteristic 
cone of Pm(/) are the vectors x G q^q]^ ^j^at a; < 0. Furthermore, the lin- 

eality space of Pm{J) is {0}. Given a polyhedron P such that lin.space P — {0}, 
it is well-known (see, e.g, [HI Chapter 8]) that any x G P can be represented 
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as X — J2^=i ^iVi + ^ where yi, . . . , y^+i are vertices of P, c G char. cone P, 
^1=1 — 1; and Ai > for all i. (That is, x is the sum of a convex combi- 
nation of some of the vertices of P and a vector in the characteristic cone of 
P.) The fact that n + 1 vertices suffice is also well-known and is a corollary to 
Caratheodory's Theorem [5] (see ^71 Chapter 7.7] for a proof of the theorem 
and the corollary). We state this result adapted to our setting as the following 
theorem. 

Theorem 5.4. Let f : R be suhmodular and let x e PAi{f )- Let N = 

n-\A\. There are vertices Hi, ... ,yN+i ofPi\i{f), coefficients Xi, X^^i e R, 
and a vector c € IR["1x-4 g^^/j ^/j^j^ 

n+l 

a; = ^ Xiyi + c, 

i=l 

c < 0, ^» = !■ firid > for each ie[N + 1]. 

We start with showing that the vertices of Pm(/) can be encoded in not too 
many bits. This is needed in the proof of Theorem 15.81 

Let m be a positive integer. Given a set of vector X C M™ we use conv(X) 
to denote the convex hull of X and cone(X) to denote 

{Xixi + ... + Xtxt \t eN,xi,...,xt e X,Xi,...,Xt>0}. 

Definition 5.5 (Facet- and vertex-complexity [H]). Let P C M™ be a polyhe- 
dron and let (p and v he positive integers. 

• P has facet-complexity at most (/) if there exists a system of linear in- 
equalities with rational coefficients that has solution set P and such that 
any equation can be encoded with at most (j) bits. If P = we require 
that (j) > m -\- 1. 

• P has vertex-complexity at most v if there exist finite sets V and E of 
rational vectors such that P = conv{V) -f cone{E) and such that each 
vector in V and E can be encoded with at most v bits. If P — % we require 
that V >m. 

Lemma 5.6 (Part of Lemma 6.2.4 in 12J). Let P C R™ be a polyhedron. If P 
has facet- complexity at most (j), then P has vertex- complexity at most Am^(f>. 

Lemma 5.7. There is a constant c such that for any suhmodular f : Ai" 7L 
the polyhedron Puif^ has vertex- complexity at most 

c ■ \A\n'^ ■ logmax(|/|). 

Proof. From the definition of -Pa/(/) it follows that PmU) has facet-complexity 
at most c- • logmax(|/|). for some constant c. The lemma now follows from 
Lemma 15.61 □ 
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Lemma [5771 tells us that the vertices of PM{f) can be encoded with not too 
many bits (that is, the size is bounded by a polynomial in n and logmax(|/|)). 
We are now ready to prove the main theorem in this section, that A4 is well- 
characterised. 

Theorem 5.8. For every fc > 3 the lattice Aik is well-characterised. 

As usual we let Ai denote an arbitrary diamond. The idea in the proof is 
that any point in PmIJ) can be represented as a convex combination of at most 
n\A\ + 1 vertices of -Pm(/) (this is Caratheodory's theorem). Furthermore, by 
Lemma 15.21 and an iterated use of Lemma 15.31 there are membership proofs for 
the vertices of Pa/(/) which can be checked in polynomial time. Hence, we get 
membership proofs for all of Pa/(/) which can be checked efficiently and by 
Theorem 14.31 we obtain the result. 

Proof. Let / : A^" — *■ Z be a submodular function and let m be some integer. 
We will show that if mintizM" /(*) — "i, then there is a proof of this fact which 
can be checked in time polynomial in n. 

We can assume that /(O^vfn) = as i i— > f{t) — fiO^n.) is submodular. 
Let N = n ■ \A\. The proof consists of a tuple m e iV + 1 vectors 

xx,. . . ,xn+i e M["1^^, for each i e [N + 1] a, sequence tj,. . . e A^" of 
tuples, and finally an integer- valued vector c e R["1^"^. To verify the proof we 
first find A = (Ai, . . . , Ajv+i) G R^+i and y G R^xA g^^j^ ^j^^^^ 

JV+l JV+l 

y<0, ^A,a;i +y = c, A > 0, and ^ A, = 1. (9) 

i=l i=l 

This can be done in time polynomial in n. Reject the proof if there are no 
solutions to ([5]). We proceed by checking that for each i 

• OaI" - *J C t2 □ . . . C t^ri-l ^ ^2r, ^ ^^^^ 

• tj, . . . , tf^ are a^j-tight, and 

• for any j £ [2n — 1] there it at most one coordinate I e [n] such that 
4(1) ^Om andt^'+'(0 = I.M. 

Reject the proof if any of these checks fail. We now want to verify that Xi G 
-Pm(/), this can be done by using Lemma [5T51 repeatedlv. For j = 1, 2, . . . , 2n — 1 
we use the algorithm in Lemma [5731 with a ~ t\ and b ~ t^"*""*"- If all invocations 
of the algorithm succeeds we can conclude that Xi G Pi\i{f), otherwise the proof 
is rejected. Finally, compute 

N+l 

c = ^ XiXi + y 
1=1 

and accept the proof if c < 0, c is unified, and c(1^ti) = /(m) = m. 
We now prove that this proof system is sound and complete. 
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Completeness (That is, if m = min^g^n f{y) then there is a proof which 
the verifier accept.) By Theorem 14.31 there is a unified integer- valued vector 
c such that c G -Pm(/), c < and m = c(ljv4™). By Theorem 15.41 there are 
vectors cci, . . . , a^AT+i such that for each i £ [N + 1] Xi is a. vertex of Puif) 
and c is the sum of a convex combination of Xx, . . . ,a;jv-|-i with coefficients 
A — (Ai, . . . , Aiv+i) and some vector y G char. cone PM{f), hence A,y is a 
solution to dH). 

As for each z e [iV + 1] the vector Xi is a vertex of PmU) it follows from 
Lemma lS?^ fand the observation that Ojv<-n, and Ijw". are ajj-tight) that there is a 
sequence of tuples t?, . . . , t?" such that Ojvt~ = t| E <i E • ■ • E <i " = 
and for each j G [2?! — 1] there is at most one I G [n] such that t[{l) = Otk and 
t\'^^{l) = 1m- It follows that this proof is accepted by the verifier. 
Soundness (That is, if there is a proof which the verifier accepts, then m = 
minyg^n f{y).) As the verifier accepted the proof it follows from Lemma 15.31 
that Xi G Pm{I) for each i G [A^ + 1]. As A and y is a solution to ^ it follows 
that c G Pnif) (it is a sum of a convex combination of some vectors contained 
in PKiif) and a vector in char. cone PM{f ))- From the acceptance of the verifier 
it also follows that c > 0, c is unified, and m = f{m) = c{lj^n). It now follows 
from Theorem 14.31 that m = miUy^M" f{v)- ^ 

In the proof system above, instead of letting the verifier solve ^ we could 
have required that A and y are given in the proof. However, it is not obvious that 
A and y can be encoded in 0(n'''+'^) bits (for some constant c). This follows from 
the approach taken above by the fact that there are polynomial-time algorithms 
for finding solutions to systems of linear inequalities and Lemma 15.71 

Note that the vectors given in the proof do not need to be vertices of Pm (/) . 
However, by using the tight tuples and by repeatedly using Lemma 15.31 we 
can verify that the given vectors are in fact contained in Pa/(/) anyway. By 
Lemma 15.21 vectors and tight tuples always exist which satisfies the conditions 
above (namely, if we chose some appropriate vertices of Pm(/))- 

The following lemma, which uses Lemma 15.31 essentially as we use it in 
Theorem 15.81 will be useful to us in Section [S] 

Lemma 5.9. Let k he some fixed positive integer. Let f : A^" —> Z be sub- 
modular and let x be a vector in Pm(/). Let ti, . . . , tm G be x-tight tuples 
such that 0_M.n — ti n . . . \Z tm ~ 1a4" a'^'^ for each j E [m — 1] there is at 
most k distinct ii,i2, ■ ■ ■ ,ik G [n] such that tj(«i) = tj{i2) = . .. = tj{ik) = Om 
and tj_|_i(ii) = tj+i(i2) = tj+i(ifc) = Im- 

For i G [m] let Ei C /(ij) such that e Cz Ei if and only if {e,x) — f{ti). 
Given c G ^.^ ^^j, ^^ ^^ n jg possible to compute max(c, y) subject 

to y £ Pnif) o,nd (e, y) = f{ti) for all i G [m] and e G Ei in time polynomial 
in n, logmax(|/|) and the encoding length of c. 

Note that we do not require that the running time depend polynomially on 

k. 

Proof. We construct a separation algorithm for the polyhedron 

[y e PM{f) I Vz G e G E, : (e, y) = f{U)} . (10) 



18 



6 FINDING THE MINIMUM VALUE 



The lemma then fohows from the equivalence of separation and optimisation 
given by the Ellipsoid algorithm. 

Given a vector y g ]R["1^^ we first test that for all i e [m] and e £ Ei we 
have (e, y) = f{ti). For each i G [m] we do this as follows: For each j S [n] such 
that ti{j) = 1_A4 the set of pairs of atoms a,b G A such that x{j,a) + x{j,b) 
is maximised must be a subset of the set of pairs of atoms a' ,b' G A such that 
y{j,a') + y{j,b') is maximised. (Otherwise there is some e £ Ei such that 
(a;,e) = f{ti) 7^ {y,G).) If this is the case then (y, e) ~ fiU) for all e G Ei if 
and only if {y, e) = f(ti) for some e £ Ei. 

Note that this test can be done in polynomial time in n and logmax(|/|) as 
TO < |A| • n. We can then use the algorithm in Lemma [Ol to test if y G PmH)- 
By combining these two tests we have a separation oracle for (fTO)) and hence 
the lemma follows. □ 



6 Finding the Minimum Value 

In this section we will show that there is an algorithm which finds the minimum 
value of a submodular / : tVI" — > Z in time polynomial in n and max(|/|). Note 
that from an algorithm which computes mintgM" /(*) one can construct an 
algorithm to find a minimiser of /, i.e., find a tuple y £ such that /(y) = 
rnint^M" f{t)- This can be done by for each x £ M minimising : A^"~^ — > K. 
defined by fx{t) = f{x,t). If min^g^n-i /^(f) = mint^M" fit), then there is 
a minimiser y G M"" to / such that y(l) = x. By iterating this procedure n 
times one finds a minimiser of /. 

We start with a high level description of the algorithm. The starting point 
is the separation problem for Pa/ (/) and the observation that G Pnif) if 
and only if mintgAi" /(*) ^ 0. Hence, given an algorithm for deciding if is 
contained in Pm(/) we can apply a binary search strategy to find a minimiser of 
/. (Note that for any c G IR the function / + c is submodular if / is submodular.) 

In each iteration i of the algorithm we maintain an upper bound Ui and 
lower bound U on mintgx-i f{t). If G Puif — [ui — li)/'2) (note that / — 
{ui - li)/2, i.e., the function /' : ^ M defined by f'{t) = f{t) - (u, - k)/2, 
is submodular if / is submodular), we iterate the algorithm with Wi+i — Ui 
and li+i = [ui — li)/2. Otherwise, if ^ PmU), we set u^+i ~ [ui — li)/2 
and lij^i = li. For an initial upper bound we can use Ui — /(Om^)- To find 
a lower bound li we can use Theorem 14.51 together with the greedy algorithm 
in Lemma [4.11 The running time of this algorithm is 0(5* • logmax(|/|) + n), 
where S is the time taken to decide if G PmU)- 

By the equivalence of separation and optimisation given by the Ellipsoid 
algorithm it is sufficient to solve the optimisation problem for Pm(/)- (The 
results we will need which are related to the Ellipsoid algorithm are given in 
Subsection 16.11 We refer the reader to [12] for an in-depth treatment of the 
theory related to this topic.) In the optimisation problem we are given c G 
qHxA g^j-^^ j^j,g supposed to solve max(c, y),y G PAf(/). To get the running 
time we are aiming for we must do this in time polynomial in n, max(|/|) and 
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the encoding length of c. 

To solve this problem our algorithm starts with a vertex of PmI/) and 
either finds an adjacent vertex with a strictly better measure or concludes that 
no such vertex exists. (This technique is called the primal-dual method, see [271 
Section 12.1].) The initial vertex is found by the greedy algorithm in Lemma HTTl 
To make this approach run in pseudo-polynomial time two parts are needed. 
The first one is that the existence of a pseudo-polynomial algorithm to go from 
vertex to a better one or conclude that no such vertex exists. We present an 
algorithm for this in Section [HUl The other part is that we must ensure that the 
algorithm makes enough progress in each iteration so that we get the bound on 
the running time we are aiming for. To this end, we prove in Section 16.31 that 
the vertices of Pm(/) are half-integral. 

This section is organised as follows: in Subsection 16.11 state some results 
we will need related to the Ellipsoid algorithm. In Subsection 16.21 we prove a 
couple of results of the structure of the vertices of Pm (/) . We also show that 
a submodular function can be turned into a strictly submodular function such 
that any minimiser of the latter is also a minimiser of the former. This will 
be useful to us in the subsequent parts of the algorithm. In Subsection 16.31 we 
prove that the vertices of Pm(/) are half-integral. Finally, in Subsection 16 .41 we 
show how we can go from one vertex of Pm(/) to a better one (if there is one) 
and how this can be used to construct an optimisation algorithm for Pm(/)- 

6.1 The Ellipsoid Algorithm 

In this subsection we present some definitions and results which are related to 
the Ellipsoid algorithm. They are all from the book [12^ . As in [12] we make the 
general assumption that for any oracle O there is an integer c such that when 
O is given input data of length n the length of the output is 0{n'^). 

Definition 6.1 (Oracle-polynomial time). An algorithm A, with access to an 
oracle O, runs in oracle-polynomial time if there is an integer c such that given 
any input of length n A makes 0(n'^) calls to O and performs 0{n'^) additional 
primitive operations. 

This is definition 6.2.2c in [12]. 

Definition 6.2 (Well-described polyhedron). A well-described polyhedron is a 
triple (P; n, (j)) where P C M" is a polyhedron with facet- complexity at most <j). 
The encoding length of (P; n, cj)) is (j) + n. 

This is definition 6.2.1 in [12]. 

Definition 6.3 (Strong optimization problem). Given a polyhedron P C R" 
and a vector c € Q", either 

• assert that P is empty, or 

• find a vector y (z P maximising (c, x) over P, or 
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• find a vector z G char, cone P such that {c, z) > 1. 

This is definition 2.1.4 from 42]. 

Definition 6.4 (Strong separation problem). Given a vector y G M", decide 
whether y ^ P , and if not, find a hyperplane that separates y from P; more 
exactly, find a vector c G R" such that (c, y) > max{(c, x) \ x d P}. 

This is a part of Theorem 6.4.9 in P] , 

Theorem 6.5. Let {P]n,(j)) be a well-described polyhedron. The strong sepa- 
ration problem and strong optimisation problem for (P; n, (p) can be solved in 
oracle-polynomial time given an oracle for the other problem. 

6.2 The Structure of the Vertices of Pa/ (/) 

The following lemma is stated in [29i Theorem 2.1] for the boolean lattice. 
Essentially the same proof works for modular lattices. We give a version of the 
lemma specialised to A^". 

Lemma 6.6. Let t,u E such that t%u and u%t, then 

p{t){2n- p{t)) + p{u){2n- p{u)) > p{tnu){2n- p{tnu)) + p{tUu){2n- p{tUu)) 

Proof. Let a = p(t □ u),f3 — p{t) — p{t n u), 7 = p{u) — p{t n u), and 6 = 
2n — p{t U u). Then the LHS is equal to 

(a + /3)(7 + 5) + (a + 7)(/3 + S) = 2a5 + 2f3j + aj + f35 + a/S + -/S 

as p is modular (i.e., p{t) + p{u) = p(t Uu) + p{t □ u)). The RHS is equal to 

a{l3 + 7 + (5) + (a + /? + 7)(5 = 2a5 + + [35 + a[3 + j5. 

Since > the lemma follows. □ 
The lemma above tells us that the function t 1-^ p{t){2n — p{t)) is strictly 
submodular. Note that if / : A^" ^ R is submodular, then / can be turned into 
a strictly submodular function /' : A^" ^ R by /'(*) = f{t) + ep{t){2n - p{t)). 
Observe that if e > is chosen small enough then any minimiser of /' is also a 
minimiser of /. Strictly submodular functions are an interesting subset of the 
submodular functions due to this observation and the following lemma. 

Lemma 6.7. Let f : A^" R be strictly submodular and let x be a vertex of 
Phiif)- Then, the x-tight tuples form a chain. 

Proof. Assume, for the sake of contradiction, that t,u € A^" are cc-tight and 
t % u and u ^ t. It follows that 

x{t) + x{u) = f{t) + f{u) > f{t n m) + /(* Uu) = x{u Uv) + x{u n v). 

(11) 
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The last equality follows from the fact that the x-tight tuples are closed under 
n and U. However, (fTTj) contradicts the supermodularity of x. □ 

Lemma 16.71 tells us that, for strictly submodular /, for each vertex x of 
PmH) the set of a;-tight tuples is a chain in As the dimension of Puif) is 
\A\n^ for every vertex x of PmH) there are \A\n linearly independent inequal- 
ities which are satisfied with equality by x. This means that for every such x 
there is a chain C t2 1^ ■ • • C *m in and linearly independent vectors 
ei, 62, . . . , e\A\n such that for each i G [\A\n] there is some e [to] such that 
Bi G I{tj(i)) and for all i € [\A\n] we have {ei,x) ~ f{tj(i))- Furthermore, x is 
the only vector which satisfies {ei, x) = f{ti) for all i G [|^|"]- 

For general (not necessarily strict) submodular functions the set of x-tight 
tuples is not necessarily a chain, but one can prove that for every vertex there is 
a chain of tuples such that some subset of the inequalities induced by the tight 
tuples characterises the vertex. That is, given the subset of inequalities induced 
by such a chain of tight tuples there is only one point in Pa//(/) which satisfies 
all the inequalities with equality. Formally we state this as the following lemma. 

Lemma 6.8. Let f : A^" M. be submodular and let x be a vertex of PM{f )- 
Then there is a chain ti C t2 C . • . □ tm in M" and linearly independent 
vectors Ci, 62, • . • , e|A|n; which are all x-tight, such that for each i G [l^ln] 
there is some j(i) G [m] such that G I(tj(^ij). 

Proof. Define f'{t) = f{t) + ep{t){2n — p{t)) and choose e > small. From 
Lemma 16.61 it follows that /' is strictly submodular. 

Let c G R["1^'^ such that x is the unique optimum to max(c, G 
PM{f)- Let x' be a vertex of Pnif') which is an optimum to max(c, y),y G 
PAiif')- From Lemma 16.71 it follows that as x' is a vertex of Pm(/') there 
are ti \Z C ... □ tm and a;'-tight ei,...,e„|^| as in the statement of 
the lemma. Let ei, . . . , e„|yi| be the rows of the matrix A and define b = 
{f{tj(i)),---J{tj{m))V and 

e = e • (p(tj-(i))(2n - . . . , p(<j(^))(2n - p{tj(^rn))))'^ ■ 

It follows that x' — A^^{b + e). We proceed by establishing two claims. 
Claim A. A-^b G Pm(/)- 

To see this assume for the sake of contradiction that A^^b ^ Puif), then 
there is some t G Ai" and e G I{t) such that eA^^b > f{t). However, 

ex' ^ eA-^{b + e) = eA-^b + eA^^e < f{t) + ep{t){2n - p{t)). (12) 

As eA~^b > fit) we can choose some (5 > such that eA^^b > f{t) + 5. By 
choosing e so that \ep{t){2n — p{t)) — eA^^e\ < (5 we get 

eA-^b + eA-^e > f{t) + ep{t){2n - p{t)) 

which contradicts We conclude that A^^b e Puif)- □ 
Claim B. (A^b,c) = {x,c). 
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Assume, for the sake of contradiction, that {A~^b,c) < {x,c). (The in- 
equahty < follows from our choice of x and Claim A.) Note that x G PmH) C 
-Pa/ (/')■ For sufficiently small e we get 

{x', c) = {A^^{b + e), c) < {x, c) 

which contradicts the optimality of a;'. □ 
We have shown that for every vertex x £ PmH) there is some vertex x' G 
-Pm(/') which satisfies some inequalities, given by the matrix A, with equality. 
As /' is strictly submodular it follows from Lemma [^771 that the a;'-tight tuples 
form a chain. By Claim A the inequalities in A also defines a point in PmH)- 
Furthermore, by Claim B this point maximises (y, c) over PM{f )- By our choice 
of c it follows that A~^b = x. The lemma follows. □ 



6.3 PmW) is Half-integral 

In this subsection we will prove that if / : A^" ^ Z is submodular, then the 
vertices of -Pa/(/) arc half- integral. 

Lemma 6.9. Let f : A4" M. be submodular and let x be a vertex of PM{f)- 
For each i G [n] there are three possibilities 

L x{i, a) — x{i, b) for all a,b A; or 

2. there is exactly one atom a' ^ A such that x(i,a') > im'aaeAx{i,a); or 

3. there is exactly one atom a' (z A such that x{i,a') < maXaeAx(i,a). 

Proof. As a; is a vertex of Puif) there is a c e mI"]^'^ such that x is the unique 
optimum to {c,y),y e PM{f). As the optimum exist it follows that c > 0. 
Assume, for the sake of contradiction, that there is a coordinate i G [n] such 
that the statement of the lemma does not hold for i. Let Ai be the atoms a' € A 
which satisfies x{i, a') — maxaeA x{i, a). Similarly, let A2 be the atoms a' £ A2 
which satisfies x{i, a') — max^g^^^^ x{i, a). Finally, let A3 — A\{Ai UA2). We 
will first prove the following claim. 

Claim. There are distinct atoms b,c ^ A and a;-tight tuples <i,t2 G A^" 
such that ti{i) = b and t2{i) — c, furthermore b E As or b,c E A2. 

If a G A3 let x' = X + 6x{i, o) for some small 6 > 0. As a; is the unique 
optimum it follows that x' ^ Pnif) and hence there is an a;-tight tuple t G 
such that t{i) — a. So if jAaj > 2, then the claim holds. Similarly, if \Ai\ > 2, 
then any for any a E A\Ai we get an a;-tight tuple t such that t{i) — a. (Again 
this follows from considering the vector x' = x + 5x{i,a).) 

So lAsl < 1 and {\Ai\ = 1 or \Ai\ > \A\ - 1). If {A^l = and {\Ai\ = 1 
or \Ai\ > \A\ — 1), then the statement of the lemma holds, so we must have 
IA3I = 1. This implies that \Ai \ = 1. 

Let Ai ~ {a}. If c{i,a) > X^beAa ^^i™ ^'^^ x' = x + 6x{i,o,) — 
^^b£A2 x(*7^) for some small 6 > 0. It follows that x' ^ Pnif) and hence 
there is an a;-tight tuple t with t{i) = a. In the other case, when c{i,a) < 
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J2bGA2 x' = X — Sxii, a) + SJ2beA2 follows that there is 

an £c-tight tuple t with e A2. □ 

Let 6 and c be the atoms in the claim above and let and t2 be the a;-tight 
tuples in the claim. As / is submodular we have 

/(ti U ta) + f{ti n ta) < /(ti) + f{t2) = + x{t2). 

From this inequality and the fact that x G PmH) it follows that 

a;(<i Uta) + nta) < 

/(tiut2) + /(tint2) < 
a;(ti) +a;(t2) < 
x{tx Ut2) + x{tx nta). 

We conclude that x{t^ U ta) + x(tx n fa) = + x(t-2). However, this leads 

to a contradiction: 

a?(ti) + x{t-2) = x{i, b) + x{i, c) + x{ti[i = Om]) + x{t-2[i = Ox]) < 
x{i, b) + x{i, c) + x{{U U t2)[i = Om]) + x{ti n ta) < 
x{ti Ut2) + x{ti nia) 

So the coordinate i cannot exist. □ 
The lemma above can be strengthened if \A\ = 3, in this case only 1 and 2 
are possible. To see this, assume that A = {oi, 02, 03} and x{i, ai) — x{i, 02) > 
x{i, as). Let x' = x + (5x(i, ai) - Sx{i, a2) (or x' = x - Sxii, ai) + 02) if 
c(i,ai) < c(i,a2)). As a;' ^ Pm(/) it follows that there is some cc-tight tuple t 
with i(j) = oi (or t{i) — 02). We can then proceed as in the proof above. 

We will need the following lemma from [T3] in our proof of the half-integrality 

0fPAf(/). 

Lemma 6.10. Let A he a m x n integral matrix satisfying 

m 

Y,\Aj\<2 

1=1 

for j G [n]. Then, for every square non-singular submatrix S of A, S^^ is 
half-integral. 

By combining Lemma 16.71 Lemma 16.91 and Lemma 16.101 we are able the 
obtain the following theorem which asserts the half- integrality of Pm(/)- 

Theorem 6.11. Let f : A^" Z be submodular. For any vertex x of P]\j{f) 
and any i G [n] and a & A we have x{i,a) G {1/2 • fc | fc G Z}. 

Proof. Let a? be a vertex of Pa/(/). By Lemma [6.81 there is a chain of a;-tight 
tuples ti □ . . . C tm and linearly independent vectors ei, . . . , e\A\n such that 
for each i G [\A\n] there is some j{i) G [to] such that G I{tj(^i^). We can 
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also assume that for i < i' we have < Let E be the matrix with rows 

ei, . . . , e|^|„, then x is the unique solution to Ex ~ b, where 



b= {f J {tj(2)) {tj{\A\n))y 



Let A = {ai, 02, ... , a|A|}- By Lemma [6.91 we can assume, without loss of gen- 
erality, that and x{i,a2) = x{i,a3) = ... = a;(i,a|^|). If x{i,ai) = 03(1,02) = 
. . . = x{i, a\A\) we can identify x{i, oi), . . . , x{i, a\A\) without changing the set 
of solutions to Ex = b, in the other case when x{i,ai) > mma£Ax{i,a) or 
x{i,ai) < maxaeAx{i,a) we can identify a;(i, 02), a;(i, a|^|) without chang- 
ing the set of solutions to Ex = b. After having identified these variables 
we get a system of linear equations, E'x' — b', which has a unique solu- 
tion. Furthermore, the solution to E'x' = b' is half-integral if and only if 
Ex = b has a half- integral solution (that is, if and only if x is half- integral). 
Let X C[n\ X {1, 2} such that for each i e [n], {i, 1) e X and (i, 2) e X if and 
only if a;(i,ai) > min^g^ a;(i, a) or x{i,ai) < max^g^ a;(i, a). We can describe 
the rows, e^, 63, . . . , e\A\n ^ ^' follows 

• if x{i,ai) = x{i,a2) = ... = x{i,a\A\), then e^.(i, 1) J2aeA ej(«,a); 

• otherwise (if a;(i,ai) > mmaeAx{i,a) or x{i,ai) < max^g^ a;(i, o)), then 
e'.{i,l) = ej{i,ai) and e'.{i,2) ^ J2aeA^a^a, ^ji^)- 

As the solution to E'x' = b' and Ex — b are equal, modulo the identification 
of some of the variables, there is a subset R = {ri,r2, ■ ■ ■ ,i"\x\} Q [l^W] with 
ri < r2 < ■ ■ ■ < r|x| such that the matrix £"', with rows {e'- \ i £ R}, has an 
inverse. Furthermore, this inverse is half-integral (that is, E""^ is half-integral) 
if and only if the solution to E'x' = b' is half-integral. 

It is easy to see that the entries of E" are contained in {0, 1,2}. Furthermore, 
if c is an arbitrary column of E" , then it is of the form (0, . . . , 0, 1, . . . , 1, 2, . . . , 2)-^ 
or (0, . . . , 0, 1, . . . , 1, 0, . . . , 0)^ (in these patterns, x, . . . ,x means that x occurs 
zero or more times). It follows that for each {i,k) £ X we have 



1^1 

1=1 



^(z,fc)-e;^(z,fc)| <2. 



(13) 



Following the proof of Theorem 1 in |10| we now define 

/I \ 

-1 



u = 







V 








1 / 



We can then express the inverse of E" as {UE")^^U . By ((T3)) and Lemma fB.lOl 
it follows that [UE")^^ is half-integral and hence E"~^ is half-integral as well, 
which implies that x is half-integral. □ 
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6.4 Finding Augmentations 

Let / : ^ Z be submodular. In this section we will show that there is an 
algorithm which decides if 6 Phiif) in time polynomial in n and max(|/|). 
The strategy of the algorithm is to use the equivalence between separation 
and optimisation given by the Ellipsoid algorithm and solve the optimisation 
problem for Pm (/) instead. In the optimisation problem we are given c £ Q["1 
and are supposed to find max(c, e Pm(/)- This problem is solved by 
iterating an augmentation step in which we are in some vertex x of PmH) and 
wish to find some vertex x', adjacent to x, such that (c, x') > (c, x). 

Let c G qNx^ and assume that we want to solve max(c, y),y E -Pm(/)- Let 
T be the set of all aj-tight tuples and let E C ^teTl{t) such that e G if and 
only if there is some t G T with e G I{t) and {e,x) = f{t). Finding a vector 
y G rHx^ such that there is some (5 > which satisfies {c,x) < {c,x + Sy) 
and X + Sy € -Pm(/) or conclude that no such vector y exists is equivalent to 
solving the linear program 

max(c, z) subject to Ve G E' : (e, z) < and (c, z) < 1. (14) 

(Here z contains the variables.) The optimum of this linear program is if a; is 
optimal and 1 otherwise. The separation problem for this polyhedron reduces 
to computing 

max(e, z). 

eeE 

Define /' : TW" Z as f'{t) = (n^ + i) . /(t) + p{t){2n - p{t)). It is not hard to 
see that a minimiser of /' is also a minimiscr of /. Furthermore, by Lemma l6.6l 
/' is strictly submodular. When minimising submodular functions we can thus 
assume that the function is strictly submodular. By Lemma 16.71 if a; is a vertex 
of P/\/(/'), then T (the a;-tight tuples) is a chain. This implies that |T| < 2n. 

Lemma 6.12. /// : A^" Z, is strictly submodular and x a vertex of PmH), 
then the linear program (|14p can he solved in time polynomial in n, logmax(|/|) 
and the encoding length of c. (Assuming that T is available to the algorithm.) 

Proof. As / is strictly submodular it follows from Lemma [6.71 that |T| < 2n. 
Hence, the separation problem for can be solved in polynomial time. By 
the equivalence of separation and optimisation given by the Ellipsoid algorithm 
it follows that (|14p can be solved in time polynomial in n, logmax(|/|) and 
the encoding length of c. (Note that even though \T\ < 2n, the number of 
inequalities in E may be exponential in n. In particular the tuple can 
induce as many as ('2') inequalities.) □ 
By the algorithm in Lemma [6. 121 we can find an optimal solution z to p4)) . 
We can use this algorithm to find adjacent vertices which are better (if there 
are any). We also need to find the largest 6 > such that x + 5z € Pm(/)- We 
construct an algorithm for this in Lemma l6.14| but first we need a lemma. 
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Lemma 6.13. Let y be an optimal solution to (|14p which is a vertex such that 
(c, y) = 1. Assume that there are ti, ^2 G A^", ti □ t2, and ei G £^n/{ii), 62 S 
E n I{t2) such that (ei,y) = /(ii) and (62, y) = fi^^)- Furthermore, assume 
that there is no u £ T such that tx \Z u \Z t-^ with any e £ E Cl I{u) and 
(e,y) — f{u). Then, there are no three distinct coordinates i,j,k G [n] such 
that ti{i) — ti{j) — ti{k) — Om and t^ii) = i2(i) = *2(fc) = ^M- 

Proof. Let E' Q E he the vectors which define tight inequahties for y. As y 
is a vertex and (c, y) = 1, it follows that the polyhedron P = {2; G 
(e, z) = f{t),e G E' , e G lit)} is one dimensional. 
Let a, /3 G K be arbitrary and define y' G rWxA 

y' = y+{a + (i)xi - axj ~ PXk- 

From the non-existence of any u G Al" such that ti \Z u n t2 and e G Enl{u), 
{e,y) = f{u) it follows that (e, y') = f{t) for all e G £",6 G /(t). However, 
this means that y' £ P and as a and /3 where arbitrary it follows that P is not 
one-dimensional. This is a contradiction and the lemma follows. □ 
The following lemma is a crucial part of our pseudo-polynomial time algo- 
rithm for SFM(Al). With the algorithm in this lemma we are able to go from 
one vertex in Pa/(/) to a better one (if there is a better one). 

Lemma 6.14. Let f : Al" 1, be a strictly submodular function. Given 
c G Q["1X'4^ vertex x of Puif), and the set of x-tight tuples T , there is an 
algorithm which is polynomial in n, logmax(|/|) and the encoding length of c 
which finds a vertex y G Pnif) such that {c,y) > {c,x) or concludes that no 
such vertex exist. If y exists the set of y -tight tuples can be computed within the 
same time bound. 

Proof. If there is such a vertex y, then the value of the optimum of the linear 
program (|14p is 1. By Lemma l6.12l this optimum y' can be found in polynomial 
time. The set of tuples T' C T which are y'-tight can be found in polynomial 
time (as |T| < 2n). Furthermore, by Lemma fG.lSl the gap between two successive 
tuples in T' is not too large. It follows from Lemma lST^ that we can find a vertex 
y of Pmif) such that (c, y) > (c, x) in polynomial time. 

It remains to find the rest of the y'-tight tuples within the stated time bound. 
By Lemma 16.131 for any consecutive tuples a, 6 in T' there are at most two 
distinct coordinates i, j G [n] such that a{i) = a(j) = 0a4 and b{i) — b{j) — 1m ■ 
We will show that for every such pair a, b in T' we can find the y-tight tuples t 
which satisfies a \Z t \Z b. To do this, for each p,q £ M, we find the minimisers 
to the submodular function fp,q defined as fp,q{x) = f{x[i = p,j = q])—y{x[i = 
p, j = q]) over the set X — {x € M" | a C a; IZ b}. As / is submodular and y is 
supermodular it follows that /' is submodular. To minimise fpq over X we can 
minimise it over at most n'^\A\'^ intervals defined by {x G Af" \ a* ^ x \Z b^} 
where a ^ a* and ^ b (there are at most n\A\ choices for a* and at most 
n\A\ choices for b*). 

Note that each of these intervals is a product of the two element lattice and 
hence this minimisation can be done with the known algorithms for minimising 
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submodular set functions. We can use this method to find all minimisers of 
fp^q in the interval we are interested in. (When we have found one minimiser 
m we iteratively minimise fp^g over the sets {x € M" \ a \Z x \Z m} and 
{x G M" \ m \Z X \Z b}.) As the y-tight tuples is a chain in M" there are only 
a polynomial number of y-tight tuples and hence this step of the algorithm runs 
in polynomial time. Hence the set of all y-tight tuples can be found within the 
stated time bound. □ 
We are now finally ready to show the existence of a pseudo-polynomial time 
separation algorithm for Pm(/)- 

Theorem 6.15. Let f : Ad" Z be submodular. It is possible to decide ifO is 
contained in PAiif) or not in time polynomial in n and max(|/|). 

Proof. By the equivalence of separation and optimisation given by the Ellipsoid 
algorithm there is an algorithm which decides if is contained in Pm(/) or 
not which makes use of an optimisation oracle for Pa/(/). The number of calls 
to the optimisation oracle is bounded by a polynomial in n and logmax(|/|), 
furthermore the objective function given to the optimisation oracle is given by a 
vector c G QI"! such that the encoding length of c is bounded by a polynomial 
in n and logmax(|/|). 

To prove the lemma it is therefore sufficient to construct an algorithm such 
that given c G Z^"!^"^ (there is no loss of generality in assuming that c is 
integral, a simple scaling of c achieves this) it solves max{y,c),y G PmH) in 
time polynomial in n, max(|/|) and the size of the encoding of c. Let f'{t) = 
(n^ + 1) ■ f{t) + p{t){2n — p{t)). By Lemma 16.61 /' is strictly submodular. 
Furthermore, it is easy to see that any minimiser of /' is also a minimiser of /. 
By Lemma [6771 each vertex x of Pm(/') is "characterised" of a chain of a;-tight 
tuples. 

The algorithm consists of a number of iterations. In iteration j a current 
vertex Xj of Pa/(/') is computed together with its associated chain Cj of Xj- 
tight tuples. The initial vertex Xq and initial chain Co is computed by the 
greedy algorithm from Lemma |4. II 

In iteration j, either xj is the optimum or there is some other vertex a;j_|-i 
such that (a;j-|_i, c) > {x, c). To find such an xj-^-i or conclude that no such ver- 
tex exists we use the algorithm from Lemma 16.141 In the case when Xj^i exists 
we also get the chain Cj-i-i of a;j_(_i-tight from the algorithm in Lemma [6. 141 

By Theorem 16.111 the vertices of PmH) are half- integral. This implies that 
(a;j_|_i,c) > {xj,c) + 1/2. So the algorithm is polynomial if we can prove 
that the optimum value is not too far from the starting point Xq. That is, 
the difference between (c, a^o) and max(c, y),y G Pm{I) should be bounded by 
a polynomial in n, max(|/|) and the encoding length of c. Note that as the 
size of the encoding of c is bounded by a polynomial in n and logmax(|/|) it 
follows that maxjg[„] |c(i, a) \ is bounded by a polynomial in n and max(|/|). 
Furthermore, as Xq is obtained by the greedy algorithm it follows that for any 
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i £ [n], a G A we have — 2max(|/|) < XQ^i, a). We now obtain the inequahty 
— 2 max(|/|) • ( max \c{i,a)\]<{c,Xo}<{c,y}< 

\ie[n]MeA ) 

max(|/|) • njAj ( max |c(i,a)|). 

\ie[n],aeA / 

From this inequahty and the fact that maxjg[„] ^tg^ |c(i, a)| is bounded by a 
polynomial in n and max(|/|) it follows that the difference between (c, cco) 
and (c, y) is bounded by a polynomial in n and max(|/|). As the objective 
function increases by at least 1/2 in each iteration this implies that the number 
of iterations is bounded by a polynomial in n and max( | / [ ) . □ 
From Theorem 16.151 we now get our desired result, a pseudo-polynomial 
time algorithm for minimising submodular functions over diamonds. The proof 
of this final step was given in Section [SI 



7 Conclusions and Open Problems 

The most obvious open problem is to find a polynomial time algorithm, as 
opposed to a pseudo-polynomial time algorithm established in this paper, for 
minimising submodular functions over diamonds. One possible approach may be 
to use some kind of scaling technique see, e.g., [161 [19]. The pseudo-polynomial 
algorithm as it is presented here is very inefficient: it consists of a nested appli- 
cation of the Ellipsoid algorithm. Usually, one layer of the Ellipsoid algorithm 
is considered to be too inefficient to be used in practise. It would clearly be 
desirable to have a simpler and more efficient minimisation algorithm. 
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